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ABSTRACT 

, Using drift scan data, a new approach to determining surface brightness profiles, 

I and techniques for detecting low surface brightness signals, we fit the light profile 

CN I of the brightest cluster galaxy (BCG) in the rich cluster Abell 1651 out to 670 h^^ 

^ ' kpc. This radius is a significant fraction of the virial radius of the cluster (2 

^ ■ Mpc; Girardi et al. 1998|) , indicating that the sizes of the BCG and the cluster are 



c/3 



T^lj- ' comparable. We find that the profile is consistent with a de Vaucouleurs profile over 

. the radial range probed. We also find that the integrated light profile of the BCG in 

Abell 1651 contributes 36% of the total cluster light within 500 kpc. Including 
I all luminous components, we obtain M/Lj r^lGO h for the cluster, which would be 

overestimated by ~20% without the BCG halo. Furthermore, the relatively red color 
of the BCG at large radii suggests that recent disruption and tidal stripping of spirals 
' and dwarf ellipticals do not contribute significantly to the halo luminosity. The color 

■ and the form of the profile are consistent with a scenario in which the BCG forms 
O ■ 

■ from filamentary collapse during the epoch of cluster formation, with relatively little 
Q^i evolution in the past 5 Gyrs ( Dubinski 19981) . We remove the BCG and other detected 

o 



galaxies from the image and construct a two dimensional surface brightness map of 



H I the cluster core. Several knots of excess emission are found, but the total diffuse 



component is constrained to contribute less than 5% of the cluster light. 



Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Abell 



I 1651) — galaxies: formation — galaxies: structure — techniques: image processing 



^Department of Astronomy and Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064 

^University of California Observatories/Lick Observatory and Department of Astronomy and Astrophysics, 
University of CaHfornia at Santa Cruz, Santa Cruz, CA 95064 

■^Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tuscon, AZ 85721 

^Department of Astronomy, University of Washington, Box 351580, Seattle, WA 98195-1580 



- 2 - 



1. Introduction 

Measuring ttie surface brightness profiles of brightest cluster galaxies (BCGs) out to large 
radii is critical both for understanding the formation of these giant galaxies and for determining 
the mass-to-light ratios (M / L) of galaxy clusters. First, the form and color of the profiles yield 
information about the dynamical state and the distribution of stellar populations ( [Malamuth &: 



Richstone 1984| , pVIerritt 198^ [Schombert 1988| [Andreon, GarilU, fc Maccagni 199^ , pubinski 



1998 ), thus constraining the accretion history of these systems. Second, measurement of the light 
contributed by the extended profile of the BCG is required for an accurate determination of the 
total cluster luminosity, which is essential for deriving an unbiased M/L. 

Despite the importance of accurate BCG profile measurements, recent determinations 
disagree. For example, from a large sample of BCGs, Schombert (1986) finds that some of these 
galaxies have a 'cD halo' — an extended component, centered on the BCG, with a significantly 
different surface brightness profile than the central region of the galaxy. In contrast, Graham et al. 
(1996) find that BCGs can be fit with single Sersic (1968) profiles, with profiles that are typically 
shallower than a de Vaucouleurs (r^/^) law. Meanwhile, the few detailed analyses of diffuse light in 
individual clusters find BCGs that are well-described by r^/'^ law profiles. Uson, Bough, & Kuhn 
(1990, 1991) show that the radial profile of the brightest cluster galaxy in Abell 2029 is consistent 
with a de Vaucouleurs model (reduced x^=0.78 in R) out to r = (aft)-*^/^ =425 h^^ kpc, and 
Scheick & Kuhn (1994) conclude that the BCG in Abell 2670 has an r^/^ profile (reduced x^=0.28 
in V) out to 230 kpc. The latter result is surprising given that this giant elliptical is classified 
by both Oemler (1973) and Schombert (1986) as a cD galaxy with a pronounced envelope starting 
at r~80 h-^ kpc. 

To unambiguously constrain BCG formation and permit accurate M /L determination, we 
aim to resolve this disagreement in the form of BCG luminosity profiles. A new study is needed 
that extends the detailed analysis techniques used in Abell 2029 and 2670 to a statistical sample 
of clusters. However, attaining the required flatfielding accuracy with pointed CCD observations 
is a computationally and observationally intensive task (|Gudehus 1989| , Uson, Bough, fc Kuhn 



1990, 199l| , [Scheick fc Kuhn 1994|) . In this paper we develop a method of studying the diffuse light 
that minimizes the required telescope time and can be used to efficiently study a large sample of 
clusters. We employ very flat drift scan data, a new approach to determining surface brightness 
profiles, and techniques for detecting low surface brightness signals that, having been developed 
for finding high-redshift clusters (Dalcanton 199f:, Zaritsky et al. 1997), are also applicable to this 
problem. 

With these tools, we perform a detailed analysis of the distribution of light in the cluster 
Abell 1651 at z=0.084. We choose this cluster to be our first because it appears to be dynamically 
relaxed, and hence is a good system with which to test our method. Specifically, X-ray observations 
show a symmetric temperature profile co-centric with the brightest cluster galaxy, with a mean 
temperature Tx = 6.1 keV ( pMarkevitch et al. 1998 ). We focus on measuring the surface brightness 
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profile of the BCG, constraining the luminosity contribution of diffuse light, and assessing the 
relative contributions of various components to the total cluster luminosity. 



2. Data and Preliminary Reductions 

Drift scan data were obtained using the Las Campanas Im telescope, the Great Circle Camera 
( Zaritsky, Schectman, &: Bredthauer 199^ ), and the Tek#5 CCD with both the Gunn i filter 



(transformed to Cousins / using Landolt standards) and a wide band filter (hereafter denoted 
as W) that roughly covers the wavelength region between B and / (see Figure 1). This W filter 
is designed to maximize the incident signal while avoiding sky emission lines in the red and 
atmospheric refraction in the blue. Individual drift scans are 2048x13000 pixels with a plate scale 
of 0.7" pixel~^ and an effective exposure time of 95 seconds. The total time required for a single 
scan, including the time spent off-source, is ~10 minutes. We have three scans of the cluster core 
in I and two in W, for total exposure times of 4.75 minutes and 3.17 minutes. Conditions during 
the observations were photometric, with seeing of 1.5". 

A key property of the data is the intrinsic flatness of drift scans. The data must have residual 
flatness variation less than 0.5%, or these variations will be the dominant source of noise in our 
derived surface brightness map. In drift scans, pixel-to-pixel variation is minimized as data are 
clocked across the chip, so sensitivity variations are a concern only perpendicular to the readout 
direction (at a level ~ 2% in our raw data). Consequently, we construct a one dimensional 
flatfield, for which the Poisson noise is reduced significantly relative to a two-dimensional Hatfield. 
Flatfielding is accomplished in two stages. In both stages we use a set of 6 (W) or 7 (!) data 
scans, to construct a median averaged flatfield. Each scan is 13,000 pixels in length and the typical 
sky level in / is ~140 counts, so the associated Poisson noise is 0.03% per column. The first 
stage of flatfielding immediately follows bias subtraction and reduces sensitivity variation from 
2% to 0.2%. This level is below the noise from other sources; however, the remaining variation 
due to the presence of objects in the scans is correlated across columns. The second flatfielding 
stage is designed to remove this residual variation. To eliminate contamination from resolved 
objects, we use the segmentation image generated by SExtractor version 2.0.15 ( [Bertin fc Arnoiits 



1996) to generate a binary mask. We convolve this binary mask with a boxcar filter to mask 
all pixels within 7" of object detection regions. The second stage flatfield image is constructed 
using only the unmasked pixels. Subsequent to this final flatfielding, all scans are flat to <0.1%, 
which corresponds to ;U/=28.4 mag arcsec"^, and any residual column-column variation is not 
discernible. 

Similar to the sensitivity variation, temporal sky variability is a one dimensional problem 
with drift scan data. This variability is a smooth feature with maximum amplitude of ~10% of 
the sky level over the length of a scan. We apply a median boxcar of size 700" (1000 pixels or 855 
h^^ kpc) to a flatfielded version of each image in which all resolved objects are masked. We find 
that a filter of this size does not cause oversubtraction of the sky near the cluster core or bright 
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stars. The entire region within 350" (500 pixels) of the cluster core is also masked as a precaution. 
We also try an alternative method in which the sky is fit with a spline of order 5. Comparison of 
the two methods yields rms variations at a level of /i/ ~ 28 mag arcsec"^, and uncertainty at the 
level of fij ~ 30.5 mag arcsec"^ in fitting the profile of the brightest cluster galaxy. 

Following bias subtraction, flatfielding, and sky subtraction, we register all images. The I- 
and VF-band data are averaged to generate a single image for each band and we also add the 
data from both bands to maximize the signal available for tracing the BCG profile at large radii. 
Adding the images also reduces uncertainty from removal of the time variable sky component. 



3. The Cluster Components 

To identify and characterize all significant sources of luminosity in the core of the cluster, 
we proceed as follows. First, we model the brightest cluster galaxy to determine the form of 
its surface brightness profile (§3.1). Next, we remove the brightest cluster galaxy and all other 
detected objects from the image and analyze the distribution of light from diffuse matter and 
faint, undetected cluster galaxies (§3.2). Finally, we use this information to assess the relative 
contribution of each of these components to the total cluster luminosity (§3.3). 



3.1. The Brightest Cluster Galaxy 

We model the brightest cluster galaxy using the IRAF routine ellipse.^ SExtractor is again 
used to generate an object mask. Previous work has demonstrated that contamination from other 



sources can significantly bias the derived profile at low surface brightness levels (Uson, Bough, & 



Kuhn 1990, l99l| . Porter, Schneider, fc Hoessel 199 1| ), and so we mask all objects except the BCG 



and all pixels within 10" of object detection regions. This procedure eliminates 30% of the image 
data from further consideration by ellipse. The masking is then augmented by manual masking of 
regions near saturated stars and bright galaxies, which excludes an additional 15% of the pixels 
in the imagej^ Our approach is similar to that utilized by Vilchez- Gomez, Pello, & Sanahuja 
(1994) in their study of the diffuse light of Abell 2390. We generate best fit models for the BCG, 
with position angle and ellipticity first allowed to vary freely, and then with these parameters 
fixed to the values given by SExtractor. The resultant surface brightness profile, as a function of 
r = (a6)^/^, is unchanged whether these parameters are fixed or allowed to vary. We choose to fix 
both parameters. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

®The regions of excess surface brightness described in §3.2 are also masked at this point, and so are not responsible 
for any observed structure in the profile. 
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To measure the surface brightness profile of the BCG at large radii, an accurate determination 
of the sky level is critical. Figure 2 illustrates the effect on an intrinsic de Vaucouleurs profile 
of error in determination of the sky level. Such an error in the background leads either to the 
truncation of the profile or to the existence of an artificial envelope. As has been noted by a 
number of authors ( de Vaucouleurs &: de Vaucouleurs 1970| , Pemler 1976 , Melnick, Hoessel, & 



White 1977 , Uson, Bough, fc Kuhn 1990, 199l| ), systematic errors in the measured sky level can 



be induced by intrinsic background variation, spatial variation in detector response (e.g. flatness 
variation for CCDs), and contamination from the outer halos of other cluster members. For this 
data set we calculate that the uncertainty in our sky determination is 0.01% {cTgky ~ 29.5 mag 
arcsec"^ in /), which will impact the form of the profile at r > 200h^^ kpc. To unambiguously 
quantify the form of the profile at larger radii, we must employ a technique that is insensitive to 
level of the background light. 

Our alternate method uses the differential change in the flux. A/, between points in the 
profile. By taking a differential measurement, we obtain a quantity that is independent of a 
constant background level, but still contains all the information present in the luminosity profile. 
One difficulty exists with this approach. If we compute the differential change only for radially 
adjacent surface brightness measurements (A/j = /j — /i+i), then both the reduced ixl) 
and the error bars for the model parameters are dependent on radial sampling density. Higher 
sampling density leads to lower xt ^^'^ larger error bars on the model parameters because the flux 
difference between adjacent points decreases, while the associated uncertainty does not. 

Our preferred method, described in the Appendix, is to compute for each point the flux 
difference relative to all points at larger radii. We define 

1 ^ 

A/.^/.-^E/.. (1) 

j=i+l 

and compare this quantity with model predictions. This approach removes the dependence of 
xf, upon sampling density, permitting robust determination of both model parameters and their 
associated uncertainties. 

Using ellipse, we measure the profile of the BCG in Abell 1651 out to r=670 h^^ kpc. Beyond 
this radius, the angular extent of the elliptical annulus used to measure the profile exceeds the 
width of the image. We show in Figure 3a the fiux differential in I+W for the BCG in Abell 
1651, and our best fit de Vaucouleurs model. A best-fit Sersic (1968) model and a model with a 
cD envelope are also shown in the residual plot, Figure 36.0 The de Vaucouleurs model fit to the 
differential profile has re = (41.7 ± 0.8)h~^ kpc (68% confidence) with x^=13.7. The large Xu is 
due to radial oscillations in the light profile (see Figure 36), which is seen in both filters. This 
radial structure is not due to using a fixed position angle and eccentricity; the same oscillations 



^The properties of the cD envelope are equivalent to the envelope observed by Schombert (1988) for Abell 2670. 
It is modelled with a de Vaucouleurs profile with re=330 kpc and Ee = 0.008'Se, galaxy 
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are seen when these parameters are permitted to vary. The best-fit Sersic model formally has 
n = 4.3 lb 0.2 with = 14.3, which is slightly shallower than a de Vaucouleurs model. However, 
the oscillations seen in the residuals may contribute to this result. For example, if the three points 
beyond 300 kpc (which appear to be on a rising part of the oscillation pattern) are excluded 
from the fit, then the best-fit model has n=3.9, and so we conclude that the Sersic index of the 
galaxy profile is consistent with n=4 (de Vaucouleurs) to within our observational uncertainty. 
We also note that the cD envelope model shown in Figure 2>h yields x^=21-6. 

Figure 4 shows the I+W composite surface brightness profile of the BCG in Abell 1651, with 
the background level fixed using the results from the differential analysis, and the vertical scale 
set such that /i(l /i^^ kpc)=0 mag arcsec"^ (see also Table 1). The de Vaucouleurs, Sersic and cD 
envelope models shown are the same as in Figure 3. A fit to the stellar PSF is also overlaid for 
comparison, demonstrating that seeing has a negligible impact on the profile. Figure 5 shows the 
/- and VF-band profiles independently, and de Vaucouleurs fits, with the effective radius fixed to 
re=41.7 h^^ kpc. Error bars represent observed rms flux variations in the data and do not include 
systematic errors. For re=41.7 /i^^ kpc, independent fits in the /- and W^-bands respectively yield 
/ie(-^)=23.55 mag arcsec"^ with x^=9.4, and /ie(W^)=24.70 mag arcsec"^ with x^=14.7. 

In addition to measuring the form of the profile, we also test for the presence of a color 
gradient. While we are able to measure the profiles out to r ~500 /i^^ kpc in both bands, for 
r >100 h-^ kpc we are wary of systematic errors that may bias the resultant colors. Consequently, 
we restrict our attention to the inner 100 kpc. We find a mild gradient in the profile 
(A(VF - /)/(Alogr) = 0.25±0.08) from 15-100 h''^ kpc, with the halo at 100 h''^ kpc redder than 
the center of the BCG by 0.2 mag. This mild gradient is consistent with other studies that have 
found shallow or no color gradients in BCGs ( pVIackie 199^ , [Garilh et al. 1997|) . For comparison, 



in their study of 17 non-BCG elliptical galaxies (-22.5<M5 <-20), Franx, Illingworth, & Heckman 
(1989) find in the mean a slightly blue radial gradient (A(i?-ii)/(Alogr) ~ —0.1). 

The color gradient of the BCG is a potentially valuable probe of the evolutionary history of 
the system. Qualitatively, if significant recent accretion of cluster galaxies has occurred in the 
halo, then we should expect the halo to be blue relative to the core of the BCG, independent of 
whether the accreted systems are spirals or ellipticals. Spirals are bluer because of ongoing star 
formation, while the relative blueness of fainter ellipticals is primarily a metallicity effect ( [Larson 



1974 , Kaufmann fc Chariot 199£ , Ferreras, Chariot, &: Silk 1999| ). For reference, we estimate 



that if one ~ L^, spiral (~5% of the total luminosity of the BCG) is accreted within the past 2 
Gyrs and deposited uniformly at 50<r<100 kpc, the halo at these radii will be -^0.1 mag 
bluer after accretion. Conversely, for formation scenarios in which there is no significant, recent 
contribution to the outer halo from tidally disrupted systems, we expect either no gradient, or a 
mild red gradient if there has been subsequent star formation in the center of the galaxy. One 
such formation scenario is demonstrated by Dubinski (1998) in a simulation that also produces a 
de Vaucouleurs profile for the BCG. Unfortunately, quantitative model predictions are currently 
lacking, and so, while our observations qualitatively agree with the early formation scenario. 
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further modelling is needed for confirmation. 



3.2. The Cluster Surface Brightness Distribution 

Is there a discernible presence of intracluster light that is not associated/co-centric with the 
BCG? To investigate this issue, we first subtract from the image the BCG using our model from 
§3.1, and all other galaxies with /-band isophotal areas larger than 175 pixels using the IRAF 
package GIM2D ([Simard 1998|, [Mar lean k Simard 19981). GIM2D generates optimal bulge+disk 



fits for each object, with an exponential disk and de Vaucouleurs profile bulge. The motivation for 
this detailed approach is to remove not only the visible central regions of these galaxies, but also 
their contribution to the total cluster light at fainter surface brightness levels. 



Subsequent to the application of GIM2D, we use FOG AS ( Jarvis &: Tyson 1981 , Valdes 1993 ) 
to remove all remaining detected objects with m/ <21 (or xn.\Y < 22.8) and replace them with 
randomly drawn, local sky pixels. The interior regions of the previously removed bright galaxies 
are also replaced by locally drawn random sky pixels as a precaution against residual galactic 
light contaminating the cluster surface brightness map. At this stage, we also mask bright stars 
and the inner 35" of the BGG. Next, we convolve the cleaned image with a 10" (12.2 kpc) 
exponential kernel to generate a smoothed two dimensional map of the core region of the cluster. 
This kernel size is chosen as a compromise between sensitivity and resolution, being sufficiently 
small to resolve individual cluster galaxies, but wide enough to probe to /i >26 mag arcsec"^. 
Residual variations persist after smoothing at an approximate rms level ///=26.65 mag arcsec"^, 
and we detect fluctuations down to ~26 mag arcsec"^. Figures 6 illustrates the process used 
to generate the surface brightness map. For comparison, the bottom panel, in which the brightest 
cluster galaxy halo is not removed, is included. 

We find that within a radius of 400 Mpc of the BCG there exist three regions of excess 
brightness with peak surface brightnesses fij < 25.75 mag arcsec"^. These regions are denoted 
as A, B, and G in the middle right panel of Figure 6. An additional ring of excess brightness 
can also be seen at the bottom of the panel surrounding a masked star, but this ring is due to 
the extended point spread function of the star. Other regions of excess brightness can also be 
seen, such as the one between the BGG and the bright star, but these fall below our fj,j = 25.75 
mag arcsec"^ threshold. Of the three peaks, B is the brightest and largest, with a peak surface 
brightness fij=25.3 mag arcsec"^. The total excess luminosity from this region corresponds to 
L/=(1.4 lb 0.5) X 10^*^X0. Region B is aligned with the semimajor axis of the BGG (see Figure 6). 

To discern whether these enhancements are diffuse in origin, we look for evidence of faint 



*Note that this is the peak surface brightness after convolution. The true peak surface brightnesses could be 
significantly brighter if the sources have spatial scales much smaller than the smoothing kernel (i.e. individual faint 
galaxies). 
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galaxies coincident with the bright regions. In region B, the detected faint galaxy population 
cannot account for the excess flux. Coupled with its large angular size, this eliminates the 
possibility that the excess is due to a higher redshift group or cluster. This region is possibly 
the remnant of a tidally disrupted cluster galaxy. Similar low surface brightness features have 
been seen m Coma ( [Gregg fc West 1998| ), and can be explained as arising from tidal processes. 



Alternately, this region could arise from a significant local enhancement of galaxies fainter than 
can be detected in our data. This interpretation seems less plausible, because it would require 
an excess relative to the expected average cluster density of at least 20 galaxies fainter than our 
detection threshold, with no corresponding overdensity of galaxies brighter than the detection 
threshold. 

In contrast to region B, inspection of regions A and C reveals a group of faint galaxies 
associated with each surface brightness excess. These galaxies have m/ >21 mag arcsec"^, and so 
were not removed from the image by FOCAS prior to smoothing. In both cases the flux associated 
with the fluctuation can be accounted for by light from these faint galaxies. Such groups could 
either be clumped dwarf galaxies associated with the cluster, or more distant clusters or groups. 
Based on the work of Gonzalez et al. (2000) detecting distant clusters, the surface density of 
background clusters is such that ~1 would be detected in this fleld, and so is consistent with 
the latter explanation. To constrain the origin of fluctuations A and C, we measure the W — I 
color of the associated faint galaxies. These galaxies are shown in Figure 7 and compared to the 
color-magnitude relation for the cluster, derived using galaxies with <20 (dashed line). The 
weighted mean color of the ensemble is — 7=0.97 (^0.3 mag bluer than the core of the BCG), 
which corresponds to i? — /~0.58. This color is consistent with the galaxies being dwarf ellipticals 
in the cluster, and inconsistent with their being giant ellipticals at higher redshift. Based on color 
alone, however, we cannot rule out the possibility that these galaxies are clustered spirals at higher 
redshift. If these are indeed cluster dwarfs, it is unclear why they are tightly clustered apart from 
any bright cluster galaxies. 



3.3. Relative Luminosity Contributions 

In §3.2 we found that a de Vaucouleurs proflle with a total luminosity of Lj = 1.17 x lO^^/i^^ 
Lq is a good fit to the BCG light profilej^ To assess the fractional contribution of the BCG to 
the total cluster light, we examine the region within r=500 kpc of the center of the BCG. 
As mentioned earlier, the location of the BCG is consistent with the kinematic center of the 
cluster as defined by X-ray data. Roughly 98% of the luminosity of the central galaxy is contained 
within this radius. We calculate the total luminosity from cluster galaxies with m/ <20.5 within 
the same region. Statistical background subtraction is used to correct for the contribution of 
galaxies not associated with the cluster. We use the region of the image further than 2 /i^^ Mpc 



^The luminosity is corrected for a galactic extinction of 0.05 magnitudes ( 



schlegel, Finkbeiner, & Davis 1996). 
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from the cluster core, covering to an area of 0.58 sq. degrees, to compute the average off-cluster 

contribution. The integrated flux for galaxies with mj <20.5 is computed in both cluster and 
off-cluster regions, with background subtraction resulting in a 25% reduction in the observed 
cluster flux. The summed, extinction-corrected luminosity of cluster galaxies, excluding the BCG, 
is Li = 2.1 X 10^2 Lq. 

To estimate the luminosity of cluster galaxies or intracluster light below this magnitude level, 
we use the surface brightness map (middle right in Figure 6). We cannot discern whether light in 
the surface brightness map is due to faint galaxies or diffuse intracluster light, but instead place 
an upper bound on the combined contribution by summing the total residual light within the 
same r=500 kpc elliptical region and background subtracting using an off-cluster region of 
the scan. To compute the mean off-cluster sky level, we use a region adjacent to this ellipse that 
extends 3000 pixels (35') in right ascension in both directions. The largest source of uncertainty 
in this measurement is due to large scale residual gradients in the background sky level, generated 
primarily by scattered light from off-image stars. Including a conservative estimate for this 
uncertainty, we calculate an excess flux of 1101^?° counts sec ^ within r=500 h ^ kpc, which 
corresponds to L= (O.T^Qy) x 10^^h~^ Lq, or 2^2% of the total cluster light in this region. 
Because the error bars arc large, we are not able to place a meaningful constraint on the faint end 
of the luminosity function. All faint end Schechter function slopes with a > —1.98 for m >20.5 
are consistent with the observed flux to within 3a. Our data are shallow, but our observed 
luminosity function has an upturn at the faint end in both bands and a slope \a\ ~ 1.7 — 1.9. This 
measurement suggests that any excess light can be explained as arising from faint cluster galaxies. 

Table 2 lists the total and fractional contributions of each component of the net cluster 
luminosity within 500 kpc. The central dominant elliptical is a key contributor to the cluster 
luminosity. The BCG fractional contribution of 36% is in good agreement with the results from 
Uson, Bough, & Kuhn (1990, 1991) for Abell 2029 and Scheick k Kuhn (1994) for Abell 2670, 
who find that the BCG's contribute 23% of the cluster light within 780 /i"^ kpc and 30±8% of the 
total cluster light, respectively. Furthermore, as in Abell 2029, diffuse light beyond that associated 
with the BCG is a negligible contributor to the total light of this system. 



4. The Mass-to-Light Ratio 

X-ray and optical spectroscopic data exist for Abell 1651, enabling us to determine the cluster 
mass and hence the /-band mass-to-light ratio. Markevitch et al. (1998) computed an emission 

weighted temperature of 6.1ib0.4 keV (90% confidence) for the cluster. Girardi et al. (1998) 

118 
96 

P = !^:!M1 = 1.04 ± 0.24 (2) 
kT 



compute a cluster velocity dispersion a=1006^lQ^ km s ^. Combining these two value yields 
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for /i=0.6, indicating that the gas and galaxies trace the potential with the same energy per unit 
mass, and suggesting that the cluster core is not far from equilibrium]^ To compute the mass 
enclosed within 500 kpc, we follow the method of Wu (1994), which assumes an isothermal gas 
distribution. This approximation is reasonable, as the temperature profile of Abell 1651 is quite 
flat ( pVIarkevitch et al. 199^ ). Evrard, Metzler, & Navarro (1996) demonstrated in their simulations 
that isothermal (3 models give an unbiased estimate of the true mass, but also demonstrated that 
there is ~20% scatter in the relation between the computed and true mass. With the isothermal 
assumption, the mass within a given projected radius, r^, is computed as a function of /3, T, a core 
radius Tc, and a 'physical size' R for the cluster, which defines the line-of-sight depth over which 
to integrate the enclosed mass. For a projected radius of 500 kpc the derived mass is only a 
weak function of rc and R. For example, changing R from 3 h^^ Mpc to 8 Mpc modifies the 
resulting mass by 3%, while changing rc from 20 h^^ kpc to 100 h^^ kpc alters the resulting mass 
by 2%. This uncertainty is significantly smaller than the scatter from use of the isothermal model. 
For concreteness, we adopt rc=50 kpc and R=5 Mpc. Including the scatter seen in the 
simulations, we compute 

M(rp < bOOh-^kpc) = (5.4 ± 1.4) x W^'^h'^MQ. (3) 
The total /-band mass-to-light ratio within this region then is 

M/Li{rp < 500h-^kpc) = (160 ± 45)h. (4) 



A key point to note is the dependence of M/Lj on the inclusion of the extended BCG halo. 
Using only the BCG magnitude returned by SExtractor, we would have missed 50% of the light 
from the BCG, and hence 18% of the cluster light. Consequently, if such extended halos are a 
generic feature of brightest cluster galaxies, as they appear to be from our work and that of Uson, 
Bough, k Kuhn (1990, 1991) and Scheick & Kuhn (1994), then any M/L ratio that fails to account 
for this light will also overestimate M/L by ~20%. This omission of luminosity exacerbates the 
cluster baryon problem for high models ( White et al. 1993| ). 



5. Formation of the BCG Halo 

What do the properties of the BCG tell us about the formation history of this system? We 
have found that the surface brightness profile of the BCG in Abell 1651 is, to first order, consistent 
with a de Vaucouleurs model, and have observed that the profile becomes mildly redder with 
increasing radius. The uniformity of the profile over such a large radial range argues for early 



^"Girardi et al. (1998) claim that this cluster has significant substructure, in conflict with the derived value of {3 
and detailed X-ray analysis ( Markevitch et al. 1998 ). However, this finding is based on the distribution of a relatively 
small sample of 30 cluster galaxies. 
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assembly of the extended halo.|^ Meanwhile, the color of the halo indicates recently accreted 
cluster galaxies do not contribute a significant fraction of the halo luminosity, and so also supports 
early formation. 

These results are consistent with the recent work of Dubinski (1998). In his hydrodynamic 
cluster simulation, the BCG is assembled at z~0.8 and significant accretion continues until z~0.4, 
after which there are no major mergers involving the BCG. Dubinski finds that a brightest cluster 
galaxy with an r^/^ profile out to 200 kpc is formed in this simulation via the merger of 
massive galaxies during filamentary collapse. He also notes that, as we observe in Abell 1651, 
'a cD galaxy envelope did not form in this system'. In addition, the brightest cluster galaxy in 
the simulation retains a fossil alignment with the filament and the galaxy distribution. Such an 
alignment with the galaxy distribution and X-ray gas has been observed in real systems ( pastry 



T968i [Carter k Metcalfe 1980l porter, Schneider, fc Hoessel 199^ [Allen et al. T995| , [Mulchaey fc 



Zabludoff 1998). In Abell 1651 the BCG is aligned with both the X-ray gas contours (Figure 8a) 
and, more marginally, with the distribution of confirmed cluster members obtained from NED 
(Figure 8b). The alignment of Clump B with the BCG can also be interpreted in this picture as 
the recent disruption of a galaxy infalling along the direction of the filament. 



6. Summary and Conclusions 

We perform a detailed analysis of the distribution of luminous matter in the galaxy cluster 
Abell 1651. We assess the relative luminosity contributions of the brightest cluster galaxy, the 
rest of the cluster galaxy population, and any diffuse, luminous matter in the cluster that is 
unassociated with the other two components. In the process, we develop and demonstrate an 
approach for studying the distribution of luminous matter in cluster cores that allows detailed 
analysis, but requires minimal telescope time (< 1 hour per cluster on a 40" telescope). This 
technique can be applied to a large sample of clusters to test whether there exists true variation 
in the form of BCG surface brightness profiles, and also whether intracluster light is a significant 
contributor in some systems. 

In the case of Abell 1651, we find that the brightest cluster galaxy contains 36% of the total 
cluster light within r=500 kpc of the center of the BCG and the cluster. Furthermore, the 
profile of the BCG is well approximated by a de Vaucouleurs model out to r=670 kpc, which 
is consistent with it being a D type galaxy (Matthews, Morgan, & Schmidt 1964, Morgan, Kayser 



fc White 1975 , Schombert 1987[) . Our current measurements extend further in radius than those of 



Uson, Bough, & Kuhn (1990, 1991) for the cluster Abell 2029 and Scheick & Kuhn (1994) for Abell 
2670, and we find close agreement with the results of these studies. We also detect a color gradient 



^^The presence of oscillations in the profile also provides a potentially interesting constraint, however, dynamical 
modelling is first needed to assess the timescale over which such structure can be maintainted in the cluster 
environment. 
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in the profile in the sense that the halo is redder than the core of the BCG. The properties of the 

brightest cluster galaxy are consistent with the scenario of filamentary collapse and formation of 
the BCG at high redshift, as evidenced by Dubinski's (1996) numerical simulations. 

We constrain the luminosity of any additional diffuse matter in the cluster to constitute less 
than ~5% of the total luminosity. Furthermore, the residual flux observed is consistent with 
arising from cluster galaxies fainter than the magnitude limit of our data. Thus, diffuse light 
beyond that of the BCG profile is negligible for this system. We do detect one distinct patch of 
excess light along the semimajor axis of the BCG (region B) that appears to be truly diffuse light, 
possibly from a tidally disrupted galaxy. Two additional, fainter peaks in the surface brightness 
are also seen, but are coincident with faint, clumped galaxies and have fluxes that are consistent 
with arising from these faint galaxies. The colors of these faint galaxies are consistent with the 
expected color of cluster dwarf ellipticals, and so a reasonable explanation is that these are groups 
of clumped cluster dwarfs. It is also possible that these galaxies are spirals in background clusters 
or groups. 

Finally, we compute the /-band mass-to-light ratio and find M/L/=(160ib45) h. We note 
that failure to include the luminosity of the BCG halo causes M/L to be overestimated by 18%. 
If this halo profile is typical of BCGs, any study of cluster mass-to-light ratios that ignores this 

significant contribution to the cluster light will systematically overestimate M/L. Omission of 
this light will also induce a systematic bias in measurements of Q,^ derived using the resultant 
mass-to-light ratios. 
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A. Appendix 

The standard technique for fitting surface brightness profiles is via minimization of the 
surface brightness data using the equation, 

where Xv is the reduced x^j ^ is the number of data points, M is the number of free parameters. 
Hi is the measured flux at point i, fi is the value of the functional fit at point i, and ai is the 
uncertainty in the flux. For a de Vaucouleurs model there are three free parameters, re, T,^, and 
the sky level. Sersic models have n as a fourth free parameter. This approach is adequate for 
determination of Vg and Sg; however, fluctuations in the interior region of the galaxy (particularly 
correlated fluctuations arising from small-scale structure in the profile) have an inordinate impact 
on the resultant background level. For example, in our data for Abell 1651, determination of the 
sky level in this fashion leads to truncation of the profile beyond r ^ 300 h"^ kpc, despite the fact 
that the surface brightness is observed to monotonically decrease to beyond r=600 kpc. For 
Sersic models, error in the background sky level also induces error in the derived n, and so this 
method provides little leverage on the form of the profile at large radii. 

To avoid the ambiguity that arises from uncertainty in the background level, we have opted 
in this paper for a differential approach to determining the form of surface brightness profiles. A 
quick way to gain intuition for this approach is to plot the flux difference between adjacent data 
points (A/j = fi — fi+i) as a function of radius and compare this with various models. However, 
as discussed in the text, use of adjacent points has the drawback that xt^ defined as 

(iV-l)-M-l^ V'^FT^ ' 

is dependent upon the radial sampling density, and so it is not possible to determine the 
uncertainty in the derived parameters. 

Instead, for each point in the profile we compute the average flux decrement between that 
point and all other points in the profile at larger radii. Mathematically, 

1 ^ 

^fi^fi-W—T. fr (A3) 

j=i+l 

This internal referencing allows us to compute a mean x^ value for each point in the profile, with 
larger separations receiving larger weighting in determination of this mean value (because a is 
roughly constant, while the flux decrement increases). The Xu equation for this method is: 



I N-l 
2 _ ^ \ ^ 

(iV - 1) - M - 1 ^ 



(A4) 
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with xt now independent of sampling density. Because of the increased leverage at large radii, 
this method permits a more robust determination of n than is possible with the standard surface 
brightness fitting technique, while yielding comparable values of r^, and Eg. 
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Table 1. Surface Brightness Profile Data 
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^The seeing during these observations was 1.5", and so data at smaller 
radii do not reflect the true profile of the galaxy. 

^The values given for I+W have been arbitrarily normalized such that 
/x(l kpc)=0 mag arcsec"^. 

'^At this radius the /-band flux was indistinguishable from the 
background level, and so we are only able to place a lower bound on the 
magnitude. 
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Table 2. Abell 1651 /-band Luminosity Budget 
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Figure Captions 

Fig. 1. — Comparison of the broad W filter utilized in this work to standard Johnson {BV) and 
Cousins {RI) filters. The red cutoff is designed to avoid night sky lines, while maximizing incident 
flux. 

Fig. 2. — The effect of a 0.01% error in the sky level (comparable to our observations) on an initial 
de Vaucouleurs profile. Error in sky estimation can lead to either truncation of the profile or an 
artificial excess of light. 

Fig. 3. — (a) log A/ as a function of radius, with A/ defined as in equation (1). The solid line is the 
best fitting de Vaucouleurs model (re=41.7 h^^ kpc). (b) Residual plot for the upper graph showing 
variation about the model fit. Also shown are the optimal n=4.3 Sersic model (long dashed) and a 
model with a de Vaucouleurs galaxy plus a de Vaucouleurs cD envelope with re=330 kpc and 
T,e=0.008T.e,gaiaxy (short dashcd). 

Fig. 4. — Composite I+W surface brightness profile. The lines are the same models as in Figure 
3. 

Fig. 5. — Bottom: /-band (circles) and H^-band (triangles) surface brightness profiles. Solid lines 
are best fit models as in Figure 4. Top: BCG color as a function of radius. A mild reddening is 
evident from 20-100 h^^ kpc. At larger radii, the derived color is sensitive to the background level, 
and so should be viewed with caution. 

Fig. 6.— ScqTicncc of /—band images illustrating the process used to generate the surface brightness 
map. Starting with the original image (upper left), first the BCG and brightest galaxies are 
modelled and removed (upper right). Next, all objects with /xj < 21 mag arcsec"^ are replaced 
with locally drawn sky pixels and masks are applied to the bright objects (middle left). Finally, 
this image is convolved with an exponential kernel to generate the SB map (middle right). The 
most significant feature in the surface brightness map is the region marked B, which has a peak 
surface brightness of fii = 25.3 mag arcsec"^. For comparison, a surface brightness map is also 
shown in which the halo of the beg is not modelled and removed (bottom). In all images, east is 
up and north is to the left. Each frame is 700" (855 kpc) on a side. 

Fig. 7. — Color-magnitude diagram for the faint galaxies coincident with excess surface brightness 
regions A and C. The dashed line is the red sequence for the cluster constructed from galaxies with 
mw < 20. 

Fig. 8. — (a) X-ray contours overlaid on the smoothed image of the BCG halo. As in Figure 6, 
the image is 855 h-^ kpc on a side. The X-ray image is from the ROSAT PSPC, with an exposure 
time of 7435 sec. To generate to contours shown, the image is smoothed with a Gaussian of 30 
arcsec FWHM; the resulting map has a resolution of about 40" (49 kpc). This resolution is 
determined fro the FWHM of several point sources lying just outside the displayed field. The lowest 
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contour level is 3— cr above the background and contour levels are displayed in 1—a increments. 

(b)Tlie distribution of all spcctroscopically confirmed members of the cluster that lie within a 15' 
(1100 h^^ kpc) radius of the BCG. The confirmed members preferentially represent the bright end 
of the cluster luminosity function. The solid line indicates the position angle of the BCG. The inner 
box denotes the region displayed in Figure 8a. East is up, north is to the left. 
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W Filter vs. Standard Filters 
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